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Abstract. 



In this paper we study a shallow water equation derivable using the Boussinesq approximation, 
which includes as two special cases, one equation discussed by Ablowitz et aJ [Stud. Appl. Math., 

§ ; 53 (1974) 249-315] and one by Hirota and Satsuma [J. Phys. Soc. Japan, 40 (1976) 611-612]. A 

catalogue of classical and nonclassical symmetry reductions, and a Painleve analysis, are given. Of 
particular interest are families of solutions found containing a rich variety of qualitative behaviours. 
rS I Indeed we exhibit and plot a wide variety of solutions all of which look like a two-soliton for t > 

C^ • but differ radically for i < 0. These families arise as nonclassical symmetry reduction solutions 

and solutions found using the singular manifold method. This example shows that nonclassical 
symmetries and the singular manifold method do not, in general, yield the same solution set. We 
also obtain symmetry reductions of the shallow water equation solvable in terms of solutions of the 
first, third and fifth Painleve equations. 

We give evidence that the variety of solutions found which exhibit "nonlinear superposition" is 
not an artefact of the equation being linearisable since the equation is solvable by inverse scattering. 
These solutions have important implications with regard to the numerical analysis for the shallow 
water equation we study, which would not be able to distinguish the solutions in an initial value 
problem since an exponentially small change in the initial conditions can result in completely 
different qualitative behaviours. 
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1 Introduction. 

In this paper we discuss the generahsed shallow water wave (gsww) equation 

A = U^xxt + aUxUxt + PUtUxx - Uxt - Uxx = 0, (1.1) 

where a and /3 are arbitrary, nonzero, constants. This equation, together with several variants, 
can be derived from the classical shallow water theory in the so-called Boussinesq approximation 
[1]. There are two special cases of this equation which have been discussed in the literature; (i), if 
a = (3 

Uxxxt + PUxUxt + PutUxx - Uxt - Uxx = 0, (1.2) 

which we shall call the swwi equation, and (ii) if a = 2/3 

Uxxxt + "^Pu^Uxt + PutUa:x " Ua:t - U^cx = 0, (1.3) 

which we shall call the swwii equation. These equations are often written in the nonlocal form 

(set Ux = v) 

Vxxt + l^vvt - Pv^d'^^vt -Vt-V:c = 0, (1.2*) 

where (9~^/) (x) = J^ f{y) dy, which was discussed by Hirota and Satsuma [2], and 

Vxxt + 2/3OTt - Pv^d'^vt -Vt-Va: = 0. (1.3*) 

which was discussed by Ablowitz et al [3] who showed that it is solvable by inverse scattering (see 
§4.2). Furthermore Ablowitz et al [3] remark that (1.3*) reduces to the celebrated Korteweg-de 
Vries (kdv) equation 

Ut + Uxxx + Quur, = 0, (1.4) 

which also is solvable by inverse scattering [4], in the long wave, small amplitude limit. Equation 
(1.3*) also has the desirable properties of the regularized long wave (rlw) equation [5,6] 

Vxxt + VV:c - Vt - V^ = 0, (1.5) 

sometimes called the Benjamin-Bona-Mahoney equation, in that it responds feebly to short waves. 
However, in contrast to (1.3*), the rlw equation (1.5) is thought not to be solvable by inverse 
scattering (cf., [7]). 

The GSWW equation (1.1) is discussed by Hietarinta [8] who shows that it can be expressed in 
Hirota's bilinear form [9] if and only if either (i), a = /3, when it reduces to (1.2), or (ii), a = 2/3, 
when it reduces to (1.3). Furthermore, as shown below, the GSWW equation (1.1) satisfies the 
necessary conditions of the Painleve tests due to Ablowitz et al [10,11] (see §2) and Weiss et al [12] 
(see §4.1) to be completely integrable if and only if a = /3 or a = 2/3. We show in §4.2 that the 
GSWW equation (1.1) is solvable by inverse scattering techniques in these two special cases. These 
results strongly suggest that the GSWW equation is completely integrable if and only if it has one 
of the two special forms (1.2) or (1.3). 

The SWWI equation (1.2) and swwii equation (1.3) arise as a reduction of several higher- 
dimensional partial differential equations which have been discussed in the literature. The swwi 
equation (1.2) arises as a reduction of: 
1. The 2 -|- 1-dimensional equation 

Uyt I Vxxxy ^VxxVy 'jVx'Uxy U, (^l.DJ 

which reduces to the KDV equation (1.4) ify = x. Boiti et al [13] developed an inverse scattering 
scheme to solve the Cauchy problem for (1.6), for initial data decaying sufficiently rapidly at 
infinity; this was formulated as a nonlocal Riemann-Hilbert problem. 
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2. The 3 + 1-dimensional equation 

lly^ ~r Uxxxy ^^xx^y ^^x^xy ^xz ^5 \^' * ) 

which was introduced by Jinibo and Miwa [14] as the second equation in the so-called 
Kadonitsev-Petviashvili hierarchy of equations; however (1.7) is not completely integrable in 
the usual sense (see [15]). 

3. The 2 + l-dimensional equation 

Utt - Uxx - Uyy + UxUxt + %%t " Uxxtt " Uyytt = 0, (1.8) 

which was introduced by Yajima et al [16] as a model of ion-acoustic waves in plasmas; Kako 

and Yajima [17] have studied "soliton interactions" for (1.8). 

The SWWii equation (1.3) arises as a reduction of the 2 -|- 1-dimensional equation 

Uxt ~r y^xxxy ^'^xx^y ^^x^xy U, l-'^-'^J 

which, like (1.6), reduces to the kdv equation (1.4) ii y = x, though note that the term Uyt in 
(1.6) is replaced by Uxt in (1.9). Bogoyaviemskii [18,19] discusses the inverse scattering method of 
solution for (1.9). 

In §§2 and 3, we find first the classical Lie group of symmetries and associated reductions of 
(1.1) and then nonclassical symmetries and reductions of (1.1). The classical method for finding 
symmetry reductions of partial differential equations is the Lie group method of infinitesimal 
transformations (cf., [20,21]). Though this method is entirely algorithmic, it often involves a large 
amount of tedious algebra and auxiliary calculations which can become virtually unmanageable if 
attempted manually, and so symbolic manipulation programs have been developed, for example in 
MACSYMA, MAPLE, MATHEMATICA, MUMATH and REDUCE, to facilitate the calculations. A survey 
of the different packages presently available and a discussion of their strengths and applications is 
given by Hereman [22]. In this paper we use the macsyma program symmgrp.max [23]. 

In recent years the nonclassical method due Bluman and Cole [24] (in the sequel referred to 
as the nonclassical method, see §3 for further details), which is also known as the "method of 
conditional symmetries" [25] or the "method of partial symmetries of the first type" [26], and 
the direct method of Clarkson and Kruskal [27] have been used to generate many new symmetry 
reductions and exact solutions for several physically significant partial differential equations that are 
not obtainable using the classical Lie method, which represents important progress (see for example 
[28,29] and references therein). Since solutions of partial differential equations asymptotically tend 
to solutions of lower-dimensional equations obtained by symmetry reduction, some of these special 
solutions will illustrate important physical phenomena. In particular, exact solutions arising from 
symmetry methods can often be effectively used to study properties such as asymptotics and "blow- 
up" (cf. [30,31]). Furthermore, explicit solutions (such as those found by symmetry methods) can 
play an important role in the design and testing of numerical integrators; these solutions provide 
an important practical check on the accuracy and reliability of such integrators (cf. [32,33]). 

There is much current interest in the determination of symmetry reductions of partial differential 
equations which reduce the equations to ordinary differential equations. Often one then checks if 
the resulting ordinary differential equation is of Painleve type, i.e., its solutions have no movable 
singularities other than poles. It appears to be the case that whenever the ordinary differential 
equation is of Painleve type then it can be solved explicitly, leading to exact solutions to the original 
equation. Conversely, if the resulting ordinary differential equation is not of Painleve type, then 
often one is unable to solve it explicitly. 

The method used to find solutions of the determining equations for the infinitesimals in both the 
classical and nonclassical case is that of Differential Grobner Bases (dgbs), defined to be a basis B 
of the differential ideal generated by the system such that every member of the ideal pseudo-reduces 
to zero with respect to B. This method provides a systematic framework for finding integr ability 
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and compatibility conditions of an overdeterniined system of partial differential equations. It avoids 
the problems of infinite loops in reduction processes, and yields, as far as is currently possible, a 
"Triangulation" of the system from which the solution set can be derived more easily [34-37]. In a 
sense, a DGB provides the maximum amount of information possible using elementary differential 
and algebraic processes in a finite time. 

In pseudo-reduction, one is allowed to multiply the expression being reduced by differential, 
that is, non-constant, coefficients of the highest derivative terms of the reducing equations. The 
reason one must do this is that on nonlinear systems, the algorithms for calculating the differential 
analogue of a Grobner Basis will not terminate if only strict reduction is allowed. What this means 
in practice is that such coefficients are assumed to be nonzero. To obtain solutions of the system 
that evaluate to zero one of these coefficients, one needs to include it with the system from the start 
of the calculation. Such a solution is called a singular integral for the obvious reason (cf., [38]). 

The major problems with the dgb method in practice are its poor complexity and expression 
swell. However, on systems where the process can be completed within reasonable limits, by which 
is meant that the length of the expressions obtained is small enough to be meaningful, the output is 
extremely useful. Comparing the determining equations for classical symmetries and a triangulation 
for that system illustrates this point; see (2.7) and (2.8) below. For nonlinear systems, dgbs have 
been used effectively to solve the determining equations for nonclassical symmetries [37,39], using 
various strategies which address the complexity problem and which minimize the number of singular 
integral cases to be considered, i.e. which minimise the differential coefficients used in the pseudo- 
reduction processes. 

A much older method of finding a basis for the ideal of a system from which formal solutions may 
be derived, due to Janet, has been implemented for linear systems [40,41]. Also for linear systems 
(and linear differential-difference systems), the differential analogue of Buchberger's algorithm [42] 
for calculating an algebraic Grobner Basis has been implemented [43]. For orthonomic systems, 
those whose members are solvable for their leading derivative term, the Reid-Wittkopf Differential 
Algebra package [44] will calculate the Standard Form of the system, the number of arbitrary 
constants and functions a formal solution depends on (see also [40]), and the formal power series 
solution to any order [35] . This package handles equations with nontrivial coefficients of the leading 
derivative terms provided maple can solve the expression (algebraically) for the leading term. One 
can then systematically go through the singular integrals using the divpivs command. 

The triangulations of the systems of determining equations for infinitesimals arising in the 
classical and nonclassical methods in this article were all performed using the maple package 
diffgrob2 [45]. This package was written specifically to handle fully nonlinear equations. All 
calculations are strictly "polynomial", that is, there is no division. Implemented there are the 
Kolchin-Ritt algorithm, the differential analogue of Buchberger's algorithm using pseudo-reduction 
instead of reduction, and extra algorithms needed to calculate a dob (as far as possible using the 
current theory) , for those cases where the Kolchin-Ritt algorithm is not sufficient [36] . Designed to 
be used interactively as well as algorithmically, the package has proved useful for solving some fully 
nonlinear systems. As yet, however, algorithmic methods for finding the most efficient orderings, 
the best method of choosing the sequence of pairs to be cross-differentiated, for deciding when to 
integrate and read off coefficients of independent functions in one of the variables, for finding the 
best change of coordinates, and so on, are still the subject of much investigation. 

The nonclassical symmetry reductions obtained for (1.2) generate a wide variety of interesting 
exact analytical solutions of the equations which we plot (using maple) in the Figures. In §4 we 
apply the Painleve test due to Weiss et al [12] to (1.1), and then obtain another family of solutions 
of (1.2) using the singular manifold method [12,46]. We then discuss the scattering problems for 
(1.2) and (1.3) and show how the arbitrary functions in the solutions obtained would appear in a 
solution of (1.2) obtained by the inverse scattering method. Finally in §5 we discuss our results. 
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2 Classical Symmetries. 

To apply the classical method to the GSWW equation (1.1) we consider the one-parameter Lie 
group of infinitesimal transformations in (x, t, u) given by 



X = X + eS,{x,t,u) + O(e^), 
i=t + eT{x,t,u) +0(e^), 
u = u + e(p{x, t, u) + O(e^), 



(2.1a) 
(2.1b) 
(2.1c) 



where e is the group parameter. Then one requires that this transformation leaves invariant the 
set 

5a = {ii(x, t) : A = 0} , (2.2) 

of solutions of (1.1). This yields an overdetermined, linear system of equations for the infinitesimals 
^{x,t,u), T{x,t,u) and (j){x,t,u). The associated Lie algebra of infinitesimal symmetries is the set 
of vector fields of the form 



Odd 

V = S,{x, t, u)— — h r(x, t,u)— + cl){x, t, u) 



dx 



dt 



du 



(2.3) 



Having determined the infinitesimals, the symmetry variables are found by solving the characteristic 
equations 

dx dt dn , ^ 

= = ^ (2.4) 

^{x,t,u) T{x,t,u) (j){x,t,u)^ 

which is equivalent to solving the invariant surface condition 



^p = ^{x, t, u)ux + t(x, t, u)ut — 0(x, t, u) = 0. 
The set 5a is invariant under the transformation (2.1) provided that 



pr(4) V (A) 



A=0 



0, 



(2.5) 



(2.6) 



where pr'^-* v is the fourth prolongation of the vector field (2.3), which is given explicitly in terms 
of ^, r and (j) (cf., [21]). This yields the following fourteen determining equations, 



0, 



0, 



ra;=0, 
(ktu = 0, 



iu 



0, 



0, 6 = 0, 

ZU <,XX U, (Pii -\- ^;; 

(3(t)t - Tt - S,x = 0, 2(3(l)xu + acpxu - l^ixx = 0, 

PVxx ~r (pxxxu v^xu — 'J) O^Yxt ^(Pxu ~r t,xx — U) 

(pxxxt y^xx (Pxt — "5 0!9^x ~r Oyxxu sxxx '^sx — U. 



(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 
(2.7e) 



These equations were calculated using the macsyma package symingrp.max [23]. 

A triangulation or standard form [34-37] (see also [40,41]) of the determining equations (2.7) 
for classical symmetries of the GSWW equation (1.1) is the following system of eight equations, 

(2.8a) 
(2.8b) 



C. = 0, 


it = 


U, t,xx — U, 


Tu 


= 0, r, = 


a4>x- lix 


= 0, 


(54>t -Tt- (px 


= 0, 


</>« + Ca; = 



from which we easily obtain the following infinitesimals, 

^ = KiX + K2, T = g{t), (f)=-Ki 



2x t 

a (3 



9{t) 



(2.9) 
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where g{t) is an arbitrary function and ki, K2 and K3 are arbitrary constants. The associated vector 
fields are: 

_ d _ f _2x _t\ d _ d _ d ( ■._ ( ^( 9 Id 

dx \ a j3 J du^ 5x ' du^ \dt j3 du 

We remark that ^^{g) shows that (1.1) is invariant under the following variable coefficient "Galilean 
transformation" 

x = x, t = g{t), u = u+[g{t)-t]/(3, (2.10) 

i.e., iiu{x,t) is a solution of (1.1), then so is u{x,t). Solving (2.4), or equivalently solving (2.5), we 
obtain two canonical symmetry reductions. 

1 '] A -P 

Case 2.1 ki 7^ 0. In this case we set ^r-r = „, , -;— , ki = 1 and kt = ks = and obtain the 

g{t) fit) dt ' ' 

symmetry reduction 

u{x, t) = f{t)w{z) + - + ^, z = xf{t), (2.11) 

a p 



where w{z) satisfies 



d'^w .d^w , ^, dw^w „ d?w , ( dw\ 



It is straightforward to show using the algorithm of Ablowitz et ai [11] that this equation is of 
Painleve-type only if either (i), a = /3 or (ii), a. = 2/3; in the Appendix it is shown that in these 
two special cases (2.12) is solvable in terms of solutions of the third Painleve equation [47] 

d^-y 1 /dy\ 1 dv o by'^ + c d 

^ - ' ^1 --— + ay3 + ^ + -, (2.13) 



dx^ y \dx ) x dx ' " ' x y 

and the fifth Painleve equation, 

^ = I- + -^\ f^V - -— + ^^ ~ ^^' lay + -1 + ^ + ^^i^±ll, (2.14) 

dx^ \2y y — \)\dx) x dx x^ \ y] x y — l' 

with a, 6, c and d constants. Hence the Painleve Conjecture [10,11] predicts that a necessary 
condition for (1.1) to be completely integrable is that (a — f3){a — 2/3) = 0, i.e., only if (1.1) has 
one of the two special forms (1.2) or (1.3). We remark that the occurrence of the third and fifth 
Painleve equations is slightly surprising since for the Boussinesq equation 

Uxx + 3(tt )xx + Uxxxx ^ Utt, (Z.lOj 

symmetry reductions reduce the equation to the first, second and fourth Painleve equations [27]. 

1 A -P 

Case 2.2 ki = 0. In this case we set — -— = — , ^2 = 1 and K3 = — 1//3 and obtain the symmetry 

reduction 

u{x,t)=w{z)+t/l5, z = x-f{t), (2.16) 

where w^z) satisfies 

d'^tf , „, dw d'^w 



dz'^ dz dz"^ 

Setting W = dw/dz and integrating twice yields 



+ (« + /3)-T--n:T = o. (2.17) 



dW^ ^ 



dz 



j +l{a + p)W^ = AW + B, 



(2.18) 



where A and B are constants of integration. This equation is equivalent to the Weierstrass elliptic 
function equation 

^) =V-52P-ff3, (2.19) 

where 52 and g^ are arbitrary constants [48]. 
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3 Nonclassical Symmetries. 

There have been several generahsations of the classical Lie group method for symmetry 
reductions. Bluman and Cole [24], in their study of symmetry reductions of the linear heat 
equation, proposed the so-called nonclassical method of group-invariant solutions. This method 
involves considerably more algebra and associated calculations than the classical Lie method. In 
fact, it has been suggested that for some partial differential equations, the calculation of these 
nonclassical reductions might be too difficult to do explicitly, especially if attempted manually 
since the associated determining equations are now an overdetermined, nonlinear system. For some 
equations such as the kdv equation (1.4), the nonclassical method does not yield any additional 
symmetry reductions to those obtained using the classical Lie method, while there are partial 
differential equations which do possess symmetry reductions not obtainable using the classical 
Lie group method. It should be emphasised that the associated vector fields arising from the 
nonclassical method do not form a vector space, still less a Lie algebra, since the invariant surface 
condition (2.5) depends upon the particular reduction. 

In the nonclassical method one requires only the subset of 5a given by 



5a,^ = {u{x,t) : A(u) = 0,V'(n) = 0} 



(3.1) 



where 5a is as defined in (2.2) and ^ = is the invariant surface condition (2.5), is invariant under 
the transformation (2.1). The usual method of applying the nonclassical method (e.g., as described 
in [25]), to the GSWW equation (1.1) involves applying the prolongation pr*-"^) v to the system of 
equations given by (1.1) and the invariant surface condition (2.5) and requiring that the resulting 
expressions vanish for u G 5a,j/), i-e.. 



pr(^) V (A) 



A=O,-0=O 



0, 



pi-l. ) v{'ip) 



A=0,V'=0 



0. 



(3.2) 



It is easily shown that 



pr(^) V(V^) = - {^uUa: + TuUt - (pu) 1p, 



which vanishes identically when -0 = without imposing any conditions upon ^, r and (p. However 
as shown in [39], this procedure for applying the nonclassical method can create difficulties, 
particularly when implemented in symbolic manipulation programs. These difficulties often arise 
for equations such as (1.1) which require the use of differential consequences of the invariant 
surface condition (2.5). In [39] we proposed an algorithm for calculating the determining 
equations associated with the nonclassical method which avoids many of the difficulties commonly 
encountered; we use that algorithm here. 

In the canonical case when r ^ we set r = 1. We omit the special case r = 0; in that case 
one obtains a single condition for (p with 424 summands, and even considering the subcase </>„ = 
leads to an equation more complex than the one we are studying. Eliminating ut, Uxt and u^xxt, 
in (1.1) using the invariant surface condition (2.5) yields 



A = 



e 



+ 'j'PuuUxUxx + 'J'PuxUxx + (puuuUx + 3(/), 

,4 



r^.T 



^x 



'•x 



xuuUx + 30, 
3 



ee 



lUx + 



xxu ^x 
uu^x^xx 



,u. 



+ (au 



^sxx'^xx '^^x'^xxx suuu'^x '^<,XUU^x 
^<,xu^x^xx su y^^^x^xxx \ ^U^^J <^Uxxxx 
1) [(px + (puUx - ixUx - iuUl - ^Uxx] + Uxx [P (0 - ^Ux 



1] = 0, (3.3) 



with t a parameter, which involves the infinitesimals ^ and that are to be determined. Now we 
apply the classical Lie algorithm to this equation using the fourth prolongation pr'^^-' v and eliminate 
Uxxxx using (3.3). This yields the following overdetermined, nonlinear system of equations for ^ 
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and (/) (contrast the classical case discussed in §2 above where they are linear). 

^u = 0, (3.4i) 

<Puu = 0, (3.4ii) 

{a + P){(t>u + ^.) = 0, (3.4iii) 

?</>*« + 3e'^:.x + 3^t^x - S^^xt - S^^(t>.u - itK = 0, (3.4iv) 

aC<A« - aitK + a^ft^tu - '^Pfff'xu - a^^(l>xu + P^^^xx - a^il + a^tix - aCCxt = 0, (3.4v) 

^^<t>xu4>xx + P^4>4>xx - ^t(t>xxx - ^^Cxx(t>xx - ^4>xx + a^4>l + ^S.(t>xxu(f>x - ^Cxxx<Px 

- '^iix4>x + it(t>x + C4>4>xxxu - S.(t><pxu + S.(j)xxxt - ^<t>xt = 0, (3.4vi) 

^^<t>u<Pxu - ^C^x4>xu - ^it(t>xu - ai^(t>x - ^(^4>xxu + PC<t>4>u - ^iixx4>u + H(t>xtu + [ii4>t 

+ mx(t> - P^t(f> + i'^ixxx + Qiixixx + ^itixx + 2eix - i^x - ^iixxt + Ct = 0, (3.4vii) 

'^iixx4>xu + P£,'^(t>xx + a^t(f>x + mx4>xxu - '^a^4>u4>x - Q^<t>lu - Ot^ckxt - iixixxx 



+ 2.^(/)j;„ + iixxx4>u - ^S.(l)u(f>XXU - ^C^XX " '^P^4>4>XU " ^^(f^XXtu + '^iix4>U " "i^^x " ^ VXU 



- (t4>u + ((t>tu + itix - CCxt + PCCxx4> + C 4>xxxu - aC(l)4>xu 

- itixxx - iixx + iixxxt + 'iit(t>xxu = 0. (3.4viii) 

These equations were calculated using the macsyma package symmgrp.max [23]. We then used the 
method of dgbs as outlined in [37,39] to solve this system. 

Case 3.1 a + /? 7^ 0. In this case it is straightforward to obtain the condition 

ex:.^'(a + /3)(3/3-2Q)=0. 

The case 3/? — 2a = leads to no solutions different from those obtained using the classical method. 
Subcase 3.1.1 ^x 7^ 0. This is the generic case which has the solution 



C=[kiX + K2]f(t), (f)=-Kif{t) 



2x Ko — t 

u + — ^ 



1 



where f{t) is an arbitrary function and Ki 7^ and K2 are arbitrary constants. These are equivalent 
to the infinitesimals (2.9) obtained using the classical method. 

Subcase 3.1.2 ^^ = 0. In this case it is easy to obtain the condition 

<^.xe'(/3-a)(a + /3) = 0. 

There are two subcases to consider. 

(i) a ^ P, S,x = 0. In this case the solution is 

e = fit), ^ = K^f{t) + 1/13, 

where f{t) is an arbitrary function and K3 is an arbitrary constant, which is equivalent to the 
infinitesimals (2.9) obtained using the classical method in the case when ki = 0. 
(ii) a = (3, ^x = ^- In this case, we obtain the following dgb for ^, cf) 

en = 0, ix=^, 0n = O, (3.5) 

i(t)xxxx -{i + l)<Pxx + PcpCpxx + Pcpl = 0, (3.6) 

(3e<Px - Pict^t + Pitd) - et = 0. (3.7) 

Thus ^ is an arbitrary function of t, and so we set ^{t) = d//di. It is easiest to integrate (3.7) first 
using the method of characteristics which yields 

^ = 2y(C)^ + i C = x + f{t). (3.8) 
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Equation (3.6) can be integrated twice with respect to x. This yields 



dt 



^^xx + \I34>'^ 



^-' 



for some arbitrary functions A(t) and ^(t). Substituting (3.8) into (3.9) yields 



(3.9) 



dt 



dC^ 



+ I3V'^ -V 



A(*)[C-/Wl + .W + i^ + ^ = o. 



We obtain an ordinary differential equation for V^(C) if we set 



X(t) = 2^4 ( ^ 



ix{t) = 2Kif{t) ( ^ 



13 dt 2/3 



+ K5, 



where K4 and K5 are arbitrary constants, yielding 

d?V 



dC 



+ PV'-V = K4C + K5- 



(3.10) 



This equation is equivalent to the first Painleve equation [47] 

— — =Qy + X, 

dx^ 



(3.11) 



if K4 7^ 0, otherwise it is equivalent to the Weierstrass elliptic function equation (2.19). Therefore 
we obtain the infinitesimals 

where C = x + f{t), f{t) is an arbitrary function and V^(C) satisfies (3.10). 

Hence solving the characteristic equations (2.4) yields the nonclassical symmetry reduction 



u{x,t)=v{C)+w{z)+t/p, C = x + f{t), z = x-f{t), 
where f{t) is an arbitrary function and v{(^) = J_ ^(Ci) dCi and w{z) satisfy 



d'^v di; d^w d^v 

dC^^ ^dCdC2~dC2 



-A, 



and 



d'^w ^^dwd'^w d'^w 

+ 2(5- — r^- ^^ = A, 



dz4 



dz dz^ dz'^ 



(3.13) 



(3.14a) 



(3.14b) 



respectively, with A a "separation" constant. Integrating (3.14) and setting V = du/dC and 
W = dw/dz, yields 

d'^V 

+ f3V^ -V = -XC + ni, (3.15a) 



and 



dC2 

d'^W 



dz2 



+ pW^-W = Xz + iJ2, 



(3.15b) 



respectively, where ni and /i2 are arbitrary constants. If A 7^ then these equations are equivalent 
to the first Painleve equation (3.11), whilst if A = then they are equivalent to the Weierstrass 
elliptic function equation (2.19). 
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In particular, if A = 0, then equations (3.15) possess the special solutions 

V{C) = -^ sech^ (kiC) + ^^^ ^ , VF(z) = ^ sech^ {k2z) + ^^^ ^ , 

where ki = ^ (1 + 4^i/3) ' and K2 = ^ (1 + 4/U2/?) • Hence we obtain the exact solution of (1.2) 
given by 

u{x, t) = —— tanh {ki [x + /(t)]} + —- tanh {k2 [x - f{t)]} 

f3 (3 

^ x{1-2kI-2kI) ^ 2f{t){4-4) ^ t ^ ^^_^g^ 

where /(t) is an arbitrary function. 

If /ii = ;U2 = then ki = K2 = -^ and (3.16) simplifies to 

uix,t) = I tanh {i [x + fit)]} + |tanh{i [x - /(t)]} + ^. (3.17) 

In Figures 1 and 2 we plot u^ with u given by (3.17) for various choices of the arbitrary function 
f{t). This is one of the simplest, nontrivial family of solutions of (1.1) with q = /?, using this 
reduction. In Figure 1, f{t) is chosen so that f{t) ~ t + to, as t ^ cx), where to is a constant. 
Consequently all the solutions plotted in Figure 1 have a similar asymptotic behaviour as t ^ 00. 
However the asymptotic behaviours as t ^ —00 are radically different. 

In the special case when /(t) = ct, then choosing ki = ^(1 + 1/c)^''^ and K2 = ^(1 — l/c)^''^ in 
(3.16) yields the two-soliton solution of (1.2) given by 




, -, X 1/2 

tanh 



, 1 \ 1/2 
c + 1 ^ 



4c 



X + ct) 



1x1/2 

+ ( tanh 



l\ 1/2 



(x — ct) 



4c 

(3.18) 
This solution is of special interest since such two-soliton solutions are normally associated with so- 
called Lie-Backlund transformations (cf., [49]) whereas (3.18) has arisen from a Lie point symmetry, 
albeit nonclassical. A plot of (3.18) for c = 2 is given in Figure 3a where it is compared to the 
so-called resonant two-soliton solution obtained using the singular manifold method in §4.1 below. 
We remark that this "decoupling" of the nonclassical symmetry reduction solution (3.13) into 
a function oi ( = x + f{t) and a function of z = x — f{t) occurs for the GSWW equation (1.1) only 
in this special case when q = /3. 

Case 3.2 a + /? = 0. Substituting (p = u9{x, t) + a{x, t) it is easy to find the condition 

GGxx — Gx = ^■ 

Thus either ^a; = or 0{x,t) = exp{a:Ai(t) + A2(t)}, where Ai(t) and A2(t) are arbitrary functions. 
In fact, unless 9x = there are no solutions. This can be shown by substituting into the equations 
the expressions 6x = OXi and 9t = (xAi^t + ^2,4)6*, to obtain, using the usual dgb techniques, 

^xx + Xl,t — S,xXl = 0. 

This can be integrated; substituting into the equations the solution, along with 9 = expjxAi (t) + 
A2(t)}, and reading off coefficients in the independent functions, the exponentials in x, 2x, and so 
on, yields various equations which lead to an inconsistency. Thus we need only consider the cases 
9 = and 9x = 0, i.e., (^„ = and (p^u = 0. 
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Subcase 3.2.1 (p^ = 0. Substituting 0„ = into the determining equations we obtain the fohowing 
DGB for the system: 

e« = o, e:. = o, 0„ = o, (3.19) 

(34>4>xx - <t>xx - P<t)l - ^(Pxxxx + ^<t>xx = 0, (3.20) 

6(1 -/?</>) +/3^'</'. + /3^0t = O. (3.21) 

Thus ^ is an arbitrary function of t and so, as in Case 3.1.2(ii), we set ,^(t) = d//dt. We integrate 
(3.21) using the method of characteristics to obtain 

(3.22) 



*->4- 


C = x-fit). 


Substituting (3.22) into (3.20) yields 




/3 


" d^r? fd7]Y 

"^dC [dcj _ 


dC4 ^ dC2 


which is not of Painleve type. Hence we obtain the infinitesimals 






-^*)-^> 



(3.23) 



(3.24) 

where C = x — f{t) and r/(C) is a solution of (3.23). These yield the (classical) symmetry reduction 
(2.16) withz = C = a;-/(t). 

Subcase 3.2.2 (j)ux = 0. In this case we obtain the solution 

i=[x + Ke]f{t), cp = -fit) [u + {2x + t + K7)/(3] + 1/(3, 

where f{t) is an arbitrary function and Kq an K7 are arbitrary constants. This is the same as the 
general case, with a = — /3. 

4 The integrability of the shallow water wave equation (1.1). 

In this section, we give first the Painleve analysis of (1.1), and use the singular manifold method 
to find another family of solutions similar, but not equivalent, to the family (3.16), for the special 
case a = p. We show that both (1.2) and (1.3) satisfy the Painleve property and are solvable 
by the inverse scattering method, suggesting that the solutions found using the nonclassical and 
singular manifold methods do not arise because there exists some transformation that linearises the 
equation (1.2). In fact, it can be seen that the arbitrary function f{t) that occurs in the families 
of solutions found arises naturally during the inverse scattering method of solution. 

4.1 The Painleve Tests. We apply the Painleve PDE test due to Weiss et al [12] to the GSWW 
equation (1.1). The Painleve Conjecture (or Painleve ODE test) as formulated by Ablowitz et al 
[10,11] asserts that every ordinary differential equation which arises as a symmetry reduction of 
a completely integrable nonlinear partial differential equation is of Painleve type, though perhaps 
only after a transformation of variables. Ablowitz et al [11] and McLeod and Olver [7] have given 
proofs of the Painleve ODE test under certain restrictions. 

Subsequently, Weiss et al [12] proposed the Painleve PDE test as a method of applying the Pain- 
leve ODE test directly to a given partial differential equation without having to consider symmetry 
reductions (which might not exist). As for the Painleve ODE test, at present there is no rigorous 
proof of the Painleve PDE test, though a partial proof can be inferred from the partial proof of 
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the Painleve ODE test due to McLeod and Olver [7]. Despite being no means foolproof the Pain- 
leve tests appear to provide a useful criterion for the identification of completely integrable partial 
differential equations. In addition to providing a valuable first test for whether a given partial 
differential equation is completely integrable, other important information can be obtained by use 
of Painleve analysis such as Backlund transformations, Lax pairs, Hirota's bilinear representation, 
special and rational solutions for completely integrable equations and special and rational solutions 
for nonintegrable equations (see, for example, [50,51] and the references therein). 

To apply the Painleve PDE test to the GSWW equation (1.1) we seek a solution in the form 



u{x,t) = yukit)cl)''+Pix,t), cP = x + ilj{t), 



fc=0 



(4.1) 



where 'ijj{t) is an arbitrary analytic function and Uk{t)^ k = 0,1,2,..., analytic functions such 
that uo ^ 0, in the neighbourhood of an arbitrary, non-characteristic movable singularity manifold 
defined by 0(x, t) = 0, and p is a constant to be determined. By leading order analysis we find that 
p = —1 and uo{t) = 12/(a + /3), provided that a + /3 / 0. In the case when a = — /3 it is routine to 
show that (1.1) is non-Painleve. We now substitute (4.1) into (1.1) to obtain from the coefficient 
of 0^-^ 

{k + l){k- l){k -A){k-6)uk = Hk{uk-i,Uk-2, ■ ■ ■ ,Uo,tlj) , 



(4.2a) 



where 



Hk = {k - 3){k - 4) 



1 + 



dV' 
'dt 



lifc_2 + (fc-4) 



du 



fc-3 



dt 



k-l 
k 



l)(j-l)[(i-2)a + (fe-j 



{k-2){k-3){k 

dip 
2)p\ujUk-j-^ 



4)- 



du 



k-l 



dt 



dt 



ik-j-l)[{j-2)a+{k-j 



2)/3]«fc--^''^' 



dt 



(4.2b) 



for k = 0, 1, 2, . . . and where we define Uk = for A; < 0. This defines Uk unless A; = 1, A; = 4 
or fc = 6 which are the so-called resonances. At each positive resonance there is a compatibility 
condition which must be identically satisfied for the expansion (4.2) to be valid, i.e., we require 
that Hi = 0, 7/4 = and Hq = for (1.1) to have a solution of the form (4.2). The compatibility 
condition Hi = is identically satisfied which implies that Ui{t) is arbitrary. Equations (4.2) with 
A; = 2 and k = 3 yield 

and 



U3 



{a - 2(5) 



'^^ 



2(a + /3)2 
respectively. The compatibility condition H/i 

I2{a - 0){a - 2(3) 
(a + /?)3 

d^fi dV' d^ip 



1 



d^V ^d^d^wi 



dt2 
yields 



/3 



dt dt2 



dV' d^tp 



+2/3- 



dt2 dt dt2 




dV' 
d^ 



d^ip 
dt^ 



0. 



Since ip is an arbitrary function this implies that 

{a - P){a - 2(3) = 0, 



(4.3) 
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is a necessary condition for (1.1) to have a solution of the form (4.2). The compatibihty condition 
Hq = is also satisfied if and only if (4.3) is satisfied. Therefore we conclude that (1.1) has a 
solution in the form (4.1) if either (i), a = (3 ox (ii), a = 2(3. These are the same conditions for 
(2.12), which arises in the classical reduction, to be of Painleve-type and for the GSWW equation 
(1.1) to expressible in Hirota's bilinear form. If (a — /3) (a — 2/3) ^ 0, then it is necessary to introduce 
a Vi{t)(j)^{x, t) In 0(x, t) term, where Vi{t) is to be determined, into the expansion (4.1) and at higher 
orders of 4>{x,t), higher and higher powers of ln(j){x,t) are required; a strong indication of non- 
Painleve behaviour. Hence the Painleve PDE test suggests that a = P and q = 2/? are the only 
integrable cases of the GSWW equation (1.1). 

Exact solutions of the swwi equation (1.2) can be obtained using the so-called singularity 
manifold method which uses truncated Painleve expansions [12,46]. If we seek a solution of (1.2) 

in the form 

6^^^(x^ 

'*^''' ^ /3 ct>{x,t) ' ^ > 

and then equate coefficients of powers of (p to zero, we find that 0(x, t) satisfies the overdetermined 
system 

(pxxxt - (pxx - 4>xt = 0, (4.5a) 

(t>t(t>xxx - ^(t>xt(t>xx - (t>l + (f>x {3(f>xxt - 4>t) = 0. (4.5b) 

(A DGB analysis of this system leads to some very complex expressions. Although it does yield 
some ordinary differential equations in x for (j) in the various subcases they appear difficult to solve.) 
Now suppose we seek a solution of these equations in the form 

(/>(x, t) = oi exp {kix + ^it} + 02 exp {k2X + ^2*} + oo, (4.6) 

where Oq, ai, 02) i^i^ f^2, f^i and H2 are constants. It is straightforward to show that equations (4.5) 
have a solution of the form (4.6) provided that fii = ki/{ki — 1), ^2 = i^2/{i^'2 ~ 1) and ki and K2 
satisfy the constraint 

k\ — K1K2 + ^2 = 3. (4.7) 

Thus we obtain the following exact solution of the swwi equation (1.2) given by 

Kit \ ( K2t 



u{x,t) 



(3 



oiKi exp <! Kix -\ — 2 7 f + '^2^2 exp <! K2a: + -^ 



Kit 1 f K2t 

oi exp <! Kix -\ — 2 7 } + 0,2 exp < K2X -\ — ^ 7 \ + ^o 



(4.8) 



provided ki and K2 satisfy (4.7). 

It should be noted that (4.8) and (3.18) are fundamentally different solutions of the swwi 
equation (1.2) as we shall now demonstrate. The general two-soliton solution of (1.2) is given by 

6 1 + Kiexp(?7i) + K2exp(?72) + Ai2(ki +K2)exp(77i +772) ., „ . 

u{x,t) = — — — , (4.9aj 

(J 1 + exp (r/i) + exp (772) + ^12 exp (r/i + 772) 

where 

kA ^ , (ki — K2)^(k? — K1K2 + to — 3) ,, ,, 

with Ki, K2, Si and 62 arbitrary constants [2] (see also [52]). The solution (4.8) is the special case 
of (4.9) with A12 = 0; Hirota and Ito [53] refer to it as being the "resonant state" where either two 
solitons fuse together after colliding with each other or a single soliton splits into two solitons (see 
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Figure 3a). On the other hand, (3.18) is the special case of (4.9) with A12 = 1 where two soUtons 
pass through each other with no phase shift as a consequence of the interaction (see Figure 3b). 
Thus whereas both solutions are asymptotically equivalent as t — > 00, they are qualitatively very 
different as t ^ —00. This shows that nonclassical method and the singular manifold method do 
not, in general, yield the same solution set. 

4.2 Inverse Scattering. The inverse scattering method, originally developed by Gardner et 
ai [4] in order to solve the kdv equation (1.4), has led to the solution of numerous physically 
significant nonlinear evolution equations, such as the nonlinear Schrodinger, Sine-Gordon, Modified 
KDV and Boussinesq equations (cf. [54]). Nonlinear evolution equations solvable by inverse 
scattering are known to possess a number of remarkable properties which appear to characterise 
the equations, including: the existence of multi-soliton solutions, an infinite number of symmetries 
and conservation laws, Backlund transformations, a Lax pair, a bi-Hamiltonian representation, a 
prolongation structure, the Hirota bilinear representation, and the Painleve property (cf. [54]). 
However, the precise relationship between these properties has yet to be rigorously established. 

There are two special cases of the GSWW equation (1.1) which have been studied from the 
inverse scattering point of view, namely the swwi equation (1.2) and the swwii equation (1.3), or 
equivalently (1.2*) and (1.3*), respectively. Hirota and Satsuma [2] studied both (1.2*) and (1.3*) 
using Hirota's bilinear technique [9]. Equation (1.3*) is known to be solvable by inverse scattering 
[3]. Several of the aforementioned properties of completely integrable equations have been derived 
for (1.2,1.3) [52,53,55-61]. 

The scattering problem for the SWWii equation (1.3) is the second order problem [3] 

i^xx + \l3u^^ = A^, (4.10) 

with associated time-dependence 

(4A - l)^Pt = (1 - Put)il)^ + \l5u^ti^, (4.11) 

where A is the constant eigenvalue, and il^xxt = V'txrr if and only if u satisfies (1.3) We note that 
(4.10) is the time-independent Schrodinger equation which is also the scattering problem for the 
KDV equation (1.4) [4]. In contrast, the scattering problem for the swwi equation (1.2) is the third 
order problem [55,59] 

i^xxx + {\(iux - 1) V'^ = Xil), (4.12) 

with associated time-dependence 

3AV't = (1 - (3ut)tpxx + Pu^ti'x- (4.13) 

We remark that (4.12) is similar to the scattering problem 

i'xxx + j{l + Qu)iP:, + j[ux- iV^d-\vt)\ V = AV' (4.14) 

which is the scattering problem for the Boussinesq equation 

Uxxxx + 3{u^)xx+Ur,a: =Utt, (4.15) 

and which has been comprehensively studied by Deift et ai [62]. 

Only the derivative Ux appears in the scattering problem (4.12) and so an arbitrary function 
of t may be added to u without affecting this. This function can be fixed by the requirement 
that ut{x,t) ^ as X ^ oo; the scattering problem (4.12) is solvable for u such that Ut{x,t) -^ 
sufficiently rapidly as x ^ oo. Moreover, the associated time-dependence, (4.13), is invariant under 
the variable-coefficient Galilean transformation (2.10). Furthermore, note that one can integrate 
the SWWI equation (1.2) once with respect to x and so introduce an arbitrary function of t. 
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5 Discussion. 

In this paper we have discussed the shallow water equation (1.1). In particular, for the special 
case of (1.1) given by the swwi equation (1.2), using the nonclassical symmetry reduction method 
originally proposed by Bluman and Cole [24], we obtained a family of solutions (3.16) which have 
a rich variety of qualitative behaviours. This is due to the freedom in the choice of the arbitrary 
function f{t). One can choose fi{t) and /2(i) such |/i(t) — /2(i)| is exponentially small as t ^ cx), 
yet fi{t) and f2{t) are quite different as t ^ — oo, so that as t ^ oo the two solutions are essentially 
the same, yet as t ^ — oo they are radically different. In Figure 1 we show that by a judicious 
choice of f{t) we can exhibit a plethora of different solutions. 

We believe that these results suggest that solving the swwi equation (1.2) numerically for 
initial conditions such as those in the solutions plotted in Figure 1 could pose some fundamental 
difficulties. An exponentially small change in the initial data yields a fundamentally different 
solution as t ^ — cxd. How can any numerical scheme in current use cope with such behaviour? 
Recently Ablowitz et al [65] have shown that the focusing nonlinear Schrodinger equation 

iut + Uxx + \u\'^u = 0, (5.1) 

exhibits numerical chaos created by small errors on the order of roundoff. The results of Ablowitz 
et al together with those given in this paper suggest that numerical analysts need to take care to 
ensure the accuracy of their programs. 

The solution (3.16) appears to be a nonlinear superposition of solutions suggesting that the 
SWWI equation (1.2) may be linearisable through a transformation to a linear partial differential 
equation, analogous to the linearisation of Burgers' equation 

Ut = Uxx + 'iuUx, (5.2) 

which is mapped to the linear heat equation through the Cole-Hopf transformation [63,64]. If 
so then the solution (3.16) could be viewed as an artefact of the fact that the swwi equation 
(1.2) is linearisable. However as illustrated in §4, the swwi equation (1.2) can be expressed as the 
compatibility condition of a third order spectral problem. Further the associated scattering problem 
(4.12) is very similar to that for the Boussinesq equation which has been thoroughly studied by 
Deift et al [62]. This strongly suggests the swwi equation (1.2) is solvable by inverse scattering. 
Additionally, as mentioned in §4, the spatial part of the inverse scattering formalism (4.12) only 
defines u up to an arbitrary additive function of t; this arbitrary function may be incorporated into 
u using the variable-coefficient Galilean transformation (2.10). 

Since the generalised shallow water equation (1.1) is invariant under the variable-coefficient 
Galilean transformation (2.10) for all a and /3, one can take any solution of the equation and using 
(2.10) generate some interesting solutions. 

Fujioka and Espinosa [69] have discussed symmetry reductions of (1.2*) using the classical Lie 
method and direct method due to Clarkson and Kruskal [27]. They claim that the classical method 
yields no symmetry reductions and that the direct method yields symmetry reductions that are a 
subset of those we obtained in §2 using the classical method. When we applied the nonclassical 
method to (1.2) in §3, we found that ^„ = 0, consequently the results of Olver [66] (see also [67,68]) 
show that the direct and nonclassical methods yield the same reductions. The difficulty Fujioka 
and Espinosa [69] appear to have experienced is with the nonlocal term in (1.2*); considering (1.2) 
rather than (1.2*) seems to be simpler. 
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Appendix 1. The solution of (2.12) 

In this appendix we show how equation (2.12), in the special case when a = (3, i.e., 
d'^w d^w ^ dw d'^w ^ d^w ^ f dw\ 

can be solved in terms of the third Painleve equation (PHI) [47] 



15 



(A.l) 



dx^ y \dx 



Idy ^ by^ + c d 

--— + ay + + -, 

X dx X y 



with a, 6, c and d constants, and also in terms of the fifth Painleve equation (PV) 



(A.2) 



dx^ 



1 1 

+ 



dy 
dx 



1 dy {y-lf / h\ cy dy{y + 1) 

X dx x^ [ y \ X y — 1 



2y y-1 
with a, 6, c and d constants. We also demonstrate how in the special case when a = 2/3, i.e. 



(A.3) 



d'^w , d^w ^ dw d^w 

z^-r+4.^^ +3/3z- 



dz^ d^s 



dz dz'^ 



+ Pw^+4p(^] =0, 



d'^w 
~d^ 



\dz j 



(A.4) 



(2.12) can be solved in terms of PIII. 

Cosgrove and Scoufis [70] consider the the equation 



X' 



d2yx2 



dA^ 



' dA / V dA ) \ dX 

dy \ dy 
A2( A— -y +^3— + A4, 

' dA / -^dA ' 



(A.5) 



where A\^ A2, A3 and A4 are constants (their equation SD-I.b (5.5)). In the general case when Ai 
and A2 are not both zero, Cosgrove and Scoufis show that (A.5) is solvable in terms of PV (A.3) 
though the transformation 



nX)^^(^^-.) -l(l-V2a)^(.-l)-^5 



iJy-1) 



+ icx^±^ + id ""'y 



where 



A = x, 



Ai = -2d, A2 = Ic^ + 2bd - d{l - V2af , 



A3 = be + ^c(l - V2af , 



A, = y 



(1 - V2af - 26 



hd 



[I - V2af + 26 



(A.6a) 
(A.6b) 

(A.7a) 
(A.7b) 



In the case when ^1 = and A2 is unrestricted, Cosgrove and Scoufis show that (A. 6) is solvable 
in terms of PIII (A.2) through the transformation 



^(^) = 1^ ("^S "^) " ^""^'^^ - H^ + 2V^)^y + ^ + 



ex dx'^ 
T&y'^ 



A = x' 



(A. 
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where 

If 21 

A2 = -jQad, ^3 = ^c (c + 2 Va) , A^ = — \ac^ - d {c + 2^) y (A.9) 

Therefore if ^i = and ^2/0 then (A.6) is solvable in terms of both PHI (A.2) and PV (A.3) 
since the special case of PV with d = can always be solved in terms of solutions of PHI [71,72]; 
Cosgrove and Scoufis [70] remark that there are infinitely many other special cases of PV and (A.6) 
that are solvable in terms of solutions of PHI, e.g., if yls = and A\ = AAiA^ in (A.6). 

To illustrate how (A.l) and (A. 4) are solvable in terms of PHI, we differentiate (A. 5) with 
respect to X to yield 



A-'Y d'Y ,fdYY AYdY ^ / dF \ ^^ As 

dX-^^dX^ = -'[dx) +XdX+^n''dX-^ +^^^+2X- 



Integrating (A.l) once yields 



d^w d'^w , /dw\ , dw 



with Bi an arbitrary constant. Now making the transformation 

^^> z Apz' ' 

and setting /3 = 9 (without loss of generality) yields (A. 10) with Ai = 0, A2 = — 16i?i/27 and 
^3 = 0. Therefore (A. 11) is solvable in terms of PV with 



3 , -3A ^ '/'' 

4 i Bi 



27^4 /-64Si^^/^ 



32^1 ' V 27 



d = 0, 



and in terms of PHI with either (i) a and b arbitrary, c = and d = 256i?i/(27a), or (ii), a and c 
arbitrary, b = — 2-y/a and d = 256Bi/{27a). 
Analogously integrating (A. 4) once yields 

^d^w d'^w dw n^ f dw\ ^( dw 1 oA „ /. x 

" dl^ + ^^ d^^ - ^ dl + t/5- ( d?) + 1" 1^'" dl - ^"' ) = ^- (^-^^l 

with B2 an arbitrary constant. Then making the transformation 

and setting /3 = 8 yields (A. 10) with ^^ = 0, 7^2 = and A^, = -\B2. Thus (A. 12) is solvable in 
terms of PHI with either (i) c and d arbitrary, o = and b = —AB2/C, or (ii) a and b arbitrary, 
c = -4^2/(6 + 2Va) and d = 0. 

We remark that (A. 5), or equations that are equivalent to (A. 5) through a Lie point 
transformation, appear in the work of Bureau [73,74], Chazy [75], Cosgrove [76,77], Jimbo [78] 
and Jimbo and Miwa [79] (see [70] for further details). 

It is well known that PHI (A.2) possesses rational solutions and one-parameter families of 
solutions expressible in terms of Bessel functions (cf., [71,80-84]) and Backlund transformations 
which map solutions of PHI into new solutions for PHI but for different values of the parameters (cf. , 
[71,72,80-85]). Starting with these known rational and one-parameter family solutions, hierarchies 
of solutions of PHI can be generated by means of the above Backlund transformations (cf., 
[86]). Analogously, it is well known that PV (A.3) possesses rational solutions and one-parameter 
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families of solutions expressible in terms of Whittaker functions (cf., [71,80-82,87,88]) and Backlund 
transformations which map solutions of the equation into new solutions with different values of the 
parameters (cf., [71,80,81,83,88-90]). Using these special exact solutions of PHI (A. 2) and PV 
(A. 3), one can construct exact solutions of (1.2) and (1.3), though we do pursue this further here. 
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Figure captions 

Figure 1. The solution (3.17) where 

(i), fit) = ¥, 

(ii), /(t) = it + exp(-t/10), 

(iii), /(t) = it+ [1 - tanht]sint, 

(iv), /(t) = |t + 2Ai(2t), 

(v), /(t) = |t + 2 exp(-tV20) sin t, 

(vi), f{t) = h + exp(-tVlOO) sin i, 

(vii), /(t) = |t + 27r-nan-i t, 

(viii), f{t) = Ut + l/t), 

(ix), /(t) = |t(l + tanht), 

where Ai(z) is the Airy function which is the solution of Ai {z) — zAi{z) = 0, satisfying 

Ai(z) ~ ivr-V2^-i/4exp (-1^3/2) as z ^ oo and Ai(z) ~ 7r-i/2|^|-i/4 cos (f |z|3/2 + lyr) 

as z ^ —CO. 

Figure 2 "Breather" solutions. The solution (3.17) where 
(i), /(t) = 2sint + 10, 
(ii), /(t) = 2sint + l, 
(iii), /(t) = 5sint + 3, 
(iv), /(t) = 5sint + l. 

Figure 3. (a) The solution (3.18) with c = 2 and (b), the solution (4.8) with k = 1.2 and 
oo = ai = 02 = 1. 



